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Figure 1: Result* using CWA. Ranges of VT are in white 
white ranges for the SR passages art in black. 



Figure 2: Results using AD. Ranges of VT are in white 
while ranges for the SR passages are in black. 



limitations. This study was done using acute leads po- 
sitioned during electrophysiologic study and has not been 
tested or validated on chronic leads. The monomorphic 
VT was reinduced during a single electrophysiology study. 
Whether changes in heart substrate or structure may cause 
changes in the intracardiac VT morphology over time has 
not been determined. 

Conclusion. We conclude that templates constructed 
from recurrent monomorphic ventricular tachycardia may 
be feasible for confirmation of the same ventricular tachy- 
cardia occurring at a later time. This may allow assessment 



of the number of times a device has delivered therapy for 
a known VT in addition to the total number of therapies 
for cardiac rhythms which may have exceeded a rate crite- 
ria and/or mophology criteria. If further study shows this 
method to be successful, it may be desirable for future anti- 
tachycardia devices to be programmed to deliver therapies ' 
based on recognition of the morphology of the VT as well a« 
its rate. In addition, VT templates may allow confirmation 
of the recurrence of known VT morphologies. 
Acknowledgement. This work was partially supported by 
NSF grant EET-8S5185 and a grant by Medtronic, Inc. 





Figure 3: Sample tracing of intracardiac electrograms during sinus rhythm (left) and ventricular tachycardia 
fright). The top channel is the atrial signal, the lower is the ventricular signal. 
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CARDIAC BIOP OTENTIAL ANALYSIS SYSTEM AND METHOD 
CROSS REFERENCE TO RELATED APPLICATION 

This application is a continuation-in-part application 
of co-pending United States Serial Ho. 107,419 filed - 
October 9, 1987, which is also assigned to the assignee of 
the present application 

i 

BACKGROUND OF THE INVENTION 

The present invention relates to a high-frequency, 
high-resolution cardiac biopotential analysis apparatus and 
method, and more particularly to a microcomputer-based 
cardiac biopotential analysis apparatus for quantitatively 
determining in a noninvasive manner; cardiac phenomena that 
can be ascertained by analyzing cardiac electrical 
activity. 

Cardiac biopotentials arise from the discharge of 
hundreds of thousands of electrically active cells, the 
signal detected at the body surface is a composite 
determined by different types of tissue, differing 
locations of that tissue, and the type of organization (or 
disorganization) of the wave front of activation. When 
transmitted to the body surface the signal is altered in 
morphology and frequency content as a result of such 
factors as body fat, ribcage size, and position of the 
heart in relation to the lungs. All these variables lead 
to challenging signal processing problems. 

Despite nearly a century of use, the conventional 
scalar electrocardiogram ("ECG") has major shortcomings. 
Its value for the diagnosis of coronary artery disease 
(CAD) is limited. It is very useful when there has been an 
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antecedent myocardial infarct (MX) ("heart attack") which 
leads to localized fibrosis extensive enough to be 
detectable. " In the absence of previous Ml, the 
abnormalities induced by CAD in the resting EC6 in an 
asymptomatic individual are of limited sensitivity and "' 
specificity. 

The scalar EC6 is of much greater value for the 
detection of active, ongoing, ischemia. Monitoring of the 
ST segment during chest pain is a reliable indicator of 
cardiac ischemia and is used diaghostically when chest pain 
spontaneously presents itself or when chest pain is 
deliberately provoked for diagnostic purposes as in 
exercise stress testing. However even for this diagnostic 
application there occur significant numbers of false 
positives often requiring £urthttr f; mora expensive 
noninvasive tests (nuclear imaging) or invasive assessment 
through the use of cardiac catheterization and coronary 
angiography. 

For the assessment of the risk of sudden cardiac death 
due to malignant ventricular arrhythmias the conventional 
ECG is of practically no use whatsoever. Twenty four hour 
continuous ambulatory monitoring of the scalar ECG ("Holter 
monitoring 11 ) is of some value in the minority of 
individuals with significant amounts of ventricular 
extrasystoles, but at considerable expense and 
inconvenience to the patient. 

Advances in computer technology have led to attempts 
to improve the diagnostic information extracted from the 
surface ECG . One such approach is the cardiointegram (CIG) 
which has been used for tlie detection of coronary diisease . 
This approach* as described in "The Cardiointegram: 
Detection of Coronary Artery Disease in Males with Chest 
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Pain and a Normal Resting Electrocardiogram", j. 
Electrocardiography. 19(3): pp. 257-267 (1986), applies a 
process of integration over the various sections of the 
H QRST signal thereby highlighting information about the 
} interrelationships of positive to negative deflections from 
the QRS to the QT interval and T wave amplitude. Using 
this technique it has been shown that coronary artery 
disease can be detected from the resting, normal ECG with a 
sensitivity and specificity slightly less than exercise 
atrmaa electrocardiography. 

More recently Abboud et. al in "High Frequency 
Electrocardiography Using an Advanced. Method of Signal 
Averaging for Non-Invasive Detection of Coronary Artery 
Disease in Patients with Normal Conventional 
Electrocardiogram", J . Electrocardioaranhy 19(4) : pp. 
371-380 (1986) showed that high frequency components of the 
ECG (150-250 Hz) averaged in the frequency domain after 
fast-Fourier transform exhibited a characteristic "zone of 
reduced amplitude" in patients with CAD. The sensitivity 
of this technique was 75%. CIG had similar sensitivity. 

A second area in which new computer based techniques 
have been applied to electrocardiography is in the 
detection of patients at risk for malignant ventricular 
arrhythmias and sudden cardiac death. Simson in "Use of 
Signals in the Terminal QRS Complex to Identify Patients 
with Ventricular Tachycardia After Myocardial Infarction", 
Clrcviatic-n.. 64(2): pp. 235-242 (1981) showed that signal 
averaging in the time domain reveals the, presence of low 
amplitude high frequency deflections in the terminal 
portion of the QRS complex, so-called "late potentials". 
These late potentials have been correlated with 
inducibility of serious arrhythmias in the 
electrophysiology (EP) lab and with an increased risk of 
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sudden death during longterm follow up of survivors of 
heart attack. However, the detection of late potentials 
has a poor^ predictive accuracy due to the problem of false 
positive tests. 

An alternative approach to the detection of risk for 
arrhythmia uses indices of the power spectrum of the signal 
averaged QRS. Cain et al in "Fast-Pouriar Transform 
Analysis of Signal -Averaged Electrocardiograms for 
Identification of Patients Prone to Sustained Ventricular 
Tachycardia", Ctrcvlrt i<?P > 69(4) : pp. 711-720 (1984) 
showed that this approach can distinguish arrhythmia 
patients from controls and correlates with inducibility of 
arrhythmias in the EP lab. Recent attempts to reproduce 
such results have met with varied success, due to 
fundamental problems in defining length of segment for PPT 
and in distinguishing the end of the QRS from noise. 
Haberl et. al in "Comparison of Frequency and Time Domain 
Analysis of Signal Averaged Electrocardiogram in Patients 
with Ventricular Tachycardia and Coronary Artery Disease: 
Methodologic Validation and Clinical Relevance. JACC . 
12(1): pp. 150-158 (1988) applied successive FPTs that are 
shifted in time to the terminal portion of the signal 
averaged QRS to address some of the problems with the Cain 
method. Neither the time-domain indices of "late 
potentials" nor the power spectral indices of Cain or 
Haberl have been shown to be influenced by drugs. This has 
limited the application of these technologies to 
diagnostics, where the problem of false positives leads to 
their use as additional procedures at additional cost . 
Judgement of therapeutic efficacy continues to require 
additional costly invasive and noninvasive procedures. 

The fundamental limitation of techniques applied to 
the ECG to dattf is tlie^ 
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electrical signal is a complex summary of spatial and 
temporal inputs and many nonlinear dynamic features should 
be expected. In particular, neural inputs to the heart 
will have significant nonlinearities . What is true in 
health is at least equally true in disease. Thus a disease 
process can be expected to lead to characteristic ' ~ 
alterations in nonlinear properties as well as linear ones. 
An ability to, quantify abnormalities in nonlinear dynamics 
would therefore be expected to enhance diagnostic power and 
improve the assessment of therapeutic efficacy. 

It is therefore a principal object of the present 
invention to provide a noninvasive system and method for 
reliably determining myocardial physiologic properties. 

Another object of the present invention is to provide 
a noninvasive system and method for quantifying linear and 
nonlinear properties of phase and energy components within 
the frequency structure of the electrocardiogram. 

A further object of the present invention is to 
provide a noninvasive system and method for diagnosing and 
quantifying coronary artery disease. 

Another object of the present invention is to provide 
a noninvasive system and method for the detection and 
quantification of myocardial ischemia in real time, for 
example as a part of intraoperative monitoring. 

Another object of the present invention is to provide 
a noninvasive system and method for the detection of 
successful reper fusion of the infarct-related artery in 
patients given thrombolytic therapy for acute myocardial 
infarction. 
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A further object of the present invention is to 
provide a noninvasive system end method for the assessment 
of coronary artery restenosis after successful percutaneous 
transluminal coronary angioplasty* 

A further object of the present invention is to * 
provide a noninvasive system and method for the 
quantification of cardiac electrical stability fixed or 
real time and the propensity for arrhythmias whether due to 
drugs, heart disease, or neural factprs. 

A further object of the present inyention is to 
provide a noninvasive system and method for the 
quantification of the extent of malignancy of cardiac 
arrhythmias. 

A further object of the present invention is to 
provide a noninvasive system and method for the 
identification of wide-complex supraventricular tachycardia 

from sustained ventricular tachycardia. 

A further object of the present invention is to 
provide a noninvasive system and method for assessing the 
efficacy of therapy for arrhythmias and sudden cardiac 
death whether that therapy is drugs or surgery. 

A further object of the present invention is to 
provide a noninvasive system and method for quantifying the 
effects of neural and humoral inputs to the heart, 
including the sympathetic and parasympathetic systems. 

A further object of the present invention is to 
provide a noninvasive system and method for evaluating pump 
function and quantifying ejection fraction. 
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It is atill another object of the present invention to 
provide a noninvasive system and method for quantifying the 
effects of ongoing organ rejection in cardiac transplant 
patients ."* 

SUMMARY OF THE INVENTION 

The cardiac biopotential analysis system and method of 
the present invention detects and quantifies the linear and 
nonlinear dynamic properties of cardiac depolarization and 
repolarization in a noninvasive manner. In so doing the 
invention provides a method for quantifying abnormalities 
in nonlinear dynamics and thereby enables without 
limitation the detection and quantification of coronary 
artery disease (CAD), myocardial ischemia, cardiac 
electrical stability, risk of malignant ventricular 
arrhythmia, site (a) of origin of malignant arrhythmias, 
extent of malignancy of arrhythmias, degree of 
antiarrhythmic drug efficacy, neural and humoral inputs to 
the heart, pump function/ ejection fraction, and ongoing 

organ rejection in cardiac transplant patients. 

• : .... 

A suitable electrode and amplifier system are used to 
acquire the cardiac electrical signal from the body surface 
of a region of interest. Very high frequency content is 
preserved by setting band pass filters at 0.05-512 Hz. 
Digital sampling is performed and digitized data is 
transmitted over a high speed serial line to a host 
microcomputer. A sinus QRST complex or a ventricular 
ectopic beat is chosen interactively as a template. Using 
standard crosscorrelation techniques a preselected number 
of complexes which match the template are extracted. 
Autobi spectral or crossbi spectral analysis is then 
performed using either an FFT approach or a parametric 
approach. 
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Zn autobispectral analysis a complex autotriple 
product array and a real autotriple product array is 
produced for a number of beats tbat match the preselected 
i> template. All of the aut obi spectral complex triple product 
arrays are then added point by point and divided by the 
total number of beats to create an average aut obi spectral 
complex triple, product array. The autobi spectral real 
triple product arrays are averaged in the same manner to 
create an average real triple product array. The magnitude 
of each averaged point in the complex autotriple product 
array is then divided by the square root of the real triple 
product array to form an autobi coherence array. An 
autobiphase array can also be produced by deriving the arc • 
tangent of the ratio of the imaginary to real part of the 
complex autotriple product array. 

In cr os sbi spectral analysis a complex cross triple 
product array and a real crosstriple product array is , 
produced for a number of successive pairs of beats that 
match the preselected template. After averaging, the 
resultant average complex crosstriple product array and 
average real crosstriple product array are used to produce 
the crossbispectrai density array, the croaSbicoherence 
array and the crossbiphase array. 

Each of the generated bi spectral arrays can contain up 
to (nf ft/2+2) *nf ft/8 data points if a nfft-point ITT is 
used. Although all, or nearly all of the 5 values at these 
points can be expected to change from normal due to 
different interventions, drugs, or disease states, in the 
preferred embodiment only those points which show the 
greatest fidelity for tracking the diagnostic determination 
in question are utiiized to create a diagnostic criterion. 
In the preferred embodiment the ensemble df points most 
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sensitive to a particular intervention or physiologic 
process is used to create one or more clinically useful 
single value indices from the computed bi spectral arrays . 
The indices are then used as diagnostic figures of merit 
for detection and quantification of coronary artery disease 
(CAD) , myocardial ischemia, cardiac electrical stability, 
risk of malignant ventricular arrhythmia, sita(s) of origin 
of malignant , arrhythmias, extent of malignancy of 
arrhythmias, degree of antiarrhythmic drug efficacy, neural 
and humoral inputs to the heart, pump function/ ejection 
fraction-, -and ongoing organ rejection in cardiac transplant 
patients. This approach makes it possible for any, even 
unskilled, operator to meaningfully interpret the output of 
the diagnostic device. 

In situations Where continuous monitbring is required, 
indices are continuously displayed on a video terminal 
enabling an operator to interactively evaluate regions of 
interest. For record keeping purposes index values and 
other pertinent variables can be sent to a hard copy output 
device or stored to magnetic storage device, such as a disk. 

These and other objects and features of the present 
invention will be more fully understood from the following 
detailed description which should be read in light of the 
accompanying drawings in which corresponding reference 
numerals refer to corresponding parts throughout the 
several views . 



BRIEF DESCRIPTION OF T HE DRAWINGS 

Fig. 1 is a schematic diagram of the components of the 
cardiac biopotential analysis system of the present 
invention; 
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rig. 2 is a schematic diagram of the 16 channel ECG 
data acquisition system, utilizing a serial interface, of 
the cardiac biopotential analysis system as shown in Fig. 1; 

Fig. 3 is a schematic diagram of the microcomputer' 
utilized by the cardiac biopotential analysis system of 
Fig. 1; 

Fig. 4 is a block diagram of the interaction of the 
various tasks performed by the system and method of the 
present invention; 

Fig. 5 is an overview flow chart of the operation of 
the system and method of the present invention; 

, Fig. 6 is a representation of the output provided by 
the system and method of the present invention; 

Fig. 7 is a flow chart of the process utilized by the 
system and method of the present invention for the 
acquisition of ZC6 data; 

Fig. 8 (a) is a diagram of a sample PQRST complex 
utilized by the system and method of the present invention; 

Figs . 8 (b) - 8 (d) are diagrams of possible extraction 
templates utilized for bispectral analysis by the method 
and system of the present invention; 

Fig. 9 is a flow chart of the frequency domain based 
steps for producing the autobi spectrum or the 
crossbi spectrum used by the system and method of the ~ 
present; 
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Fig. 10 is a flow chart of the parametric based steps 
for producing the autobi spectrum or the crossbi spectrum 
used by the system and method of the present invention; 

Fig. 11 are sample diagrams of autobi spectral arrays 
of the QRS complex generated by the system and method of 
the present invention; 

Fig. 12 is a flow chart of the steps used to generate 
diagnostic indices by the system and method of the present 
invention; , 

Fig. 13 are diagrams of several autbbicoherence arrays 
developed during the determination of the clinical 
reference arrays; 

Figs. 14(a) - 14 (b) are tables of the sample indices 
for normal subjects and for CAD subjects, respectively, 
generated by the system and method of the present invention 
in connection with the establishment of clinical reference 

arrays; 

Fig. 15 are diagrams showing sample bicoherence values 
when the system and method of the present invention is used 
to identify patients at risk bf malignant arrhythmias; 

Figs. 16(a) - 16(b) are tables of sample indices fbr 
MI subjects and VT/VF subjects respectively; 

Figs. 17(a) - 17(b) are graphs of coefficients 
generated by the system and method of the present invention. 
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DETAILED DES CRIPTION OF THE PREFERRED EMBODIMFN1TR 

Referring to Fig. 1, the apparatus of the present 
invention is connected to a patient 100 through a set .of 
surface electrodes using a standard limb, precordial and 
orthogonal placement protocol. The EC6 signals are picked 
up by the electrodes and transmitted over a patient cable 
102 to a 16 channel ECG data acquisition system 104 with a 
serial interface'. 

The data acquisition system 104 filters, amplifies and 
digitizes the ICG waveforms and sends the digitized data to 
a microcomputer 110 via a high speed synchronous serial 
line 106. In addition, the serial line 106 can be used to 
download filtering, gain and sampling rate instructions to 
the data acquisition unit 104. 

The microcomputer 110 processes the serial data stream 
in order to generate all computed data arrays. These 
arrays are then used in conjunction with predetermined 
reference arrays derived from clinical studies to produce 
diagnostic indices which indicate the status of the 
patient. These indices are displayed on the graphics 
display 108. Printed output of the diagnostic index is 
also available on the hard copy output device 116 which is 
connected to the microcomputer 110. Interaction between 
the operator and the system is provided by means of a 
keyboard 112 and pointing device 114 with feedback on the 
graphics display 10S ; : 

The 16 channel data acquisition system 110 is shown in 
greater detail i* rig. 2. The ECG surface potential, 
detected by a surface electrode mounted on the patient 100, 
plil'ttsT through " an eleictrosurgery^ protection- circuit- 2 02 , a 
defibrillator protection circuit 206 and an amplifier/ 
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filter circuit 208 bafore baing passed on to the 
multi-channel analog to digital converter 210. 

The jalectrosurgery protection circuit 202 includes a 
radio frequency (rf) filter, which limits the rf current 
through the patient leads 101 to less than 100 micro 
amperes and thus protects the patient 100 from rf burns and 
protects the amplifiers 36 from damage resulting from 
exceeding the absolute maximum input voltage specified by 
the manufacturer. This circuit can be an LC section 
circuit consisting of a generic inductor connected In 
series, to a generic capacitor which is then connected to 
ground. 

The defibrillator protection .circuit 206 limits the 
voltage to the amplifiers 208 to a safe level when a 
defibrillator is applied to the patient 100 and discharged. 
This circuit preferably includes a neon light bulb and/or a 
parallel variable resistor connected in series to a 
_ grounded resistor. 

The amplifier/ filter circuitry 208 is controlled by a 
microprocessor 216 for default gain and filtering levels or 
alternate gain and filtering levels as requested by the 
operator. Preferred gain and filtering settings are 
discussed below. This circuitry 208 includes three stages: 
the first is a pre-amplifiar stage that can be assembled 
using a wide variety of high impedance pre-amplif iers such 
as those sold by National Semiconductor, Sunnyvale 
California; the second is a programmable filters stage 
which can utilize filters sold by Frequency Devices, 
Haverhill Massachusetts; the third stage is a programmable 
amplifiers stage which can be assembled from operational 
amplifiers used in conjunction with a multiplying digital 
to analog (D/A) converter, both of which components are 
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available from National S ami conductor . The multiplying D/A 
is used to sat the gain to the appropriate levels requested 
by the microprocessor 216. 

The high impedance pre-amplif ier of each channel will 
saturate to either the positive or negative supply voltage 
if the input of the pre-amplif ier is not terminated. This 
will lead to large positive value or a large negative value 
at the output of amplifier section 208. Such values will 
be used to identify lead failure. 

The output of all 16 channels of the amplifier/ filter 
208 is fed to a multi-channel analog, to digital converter 
(A/D) 210 which is under microprocessor 216 control for 
sampling rate settings. . The analog, signals are converted 
to digital data format suitable for input to a computer. : 
A/D converters sold by Analog Devices, Norwood, 
Massachusetts can be used for this purpose. 

The multi -channel A/D converter 210 is optically 
coupled to data bus 215 by optical isolator 214. All 
control lines to the amplifiers/ filters 208 and the A/D 
convertor 210 are also optically, isolated by optical 
isolator 212. Any known optical isolators can be used for 
this purpose. 

All DC power lines going to the amplifiers/ £ liters 
208 and A/D convertor 210 are also isolated from the AC 
power line with a DC/DC convertor 204 in order to provide 
complete patient isolation from grounds- DC/DC converters 
available from Burr Brown can be used for this purpose. 

The basic instructions fop controlling operation of 
the microprocessor 216 are stored in a read only memory 
(ROM) 218v The- random- access- memory (RAM) 219 is used as a 
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buffer memory for data, and a portion of the RAM 219 can 
also be used as program memory when a control program is 
being downloaded from the microcomputer 110. 

Serial interface 220 operates under the control of the 
microprocessor 216. The serial interface 220 is optically 
coupled with optical' isolators 222 to high speed 
synchronous serial drivers 224 to provide a synchronous 
serial link 106 between the 16 channel data acquisition 
system 104 and, any compatible high speed synchronous serial 
interface card on any computer. The serial lines are 
isolated by optical isolators 222 and DC/OC convertor 204 
to provide increased patient safety and to protect the host 
computer 110' from any transients. 

The host or microcomputer 110 of Fig. 1 is shown in 
greater detail in Fig. 3. The entire microcomputer system 
runs tinder control of a microprocessor 302 with the program 
memory for the microprocessor being stored in ROM 304. The 
RAM 306 is used for storage of intermediate data. The mass 
storage device 308 is used for storing clinical databases 
as well as archiving patient data. 

In a preferred embodiment, the microcomputer 110 
contains an array processor 310 (such as the Vortex sold by 
SKY of Lowell, MA) on which complex arithmetic calculations 
can be performed on entire arrays of data simultaneously. 
The. preferred embodiment also includes a math coprocessor 
312 which is connected directly to microprocessor 302. The 
math coprocessor 312 is used for scalar and graphic 
calculations while the array processor 310 is used to 
calculate bispectral and other data vectors. 

A graphics controller 314 operating under program 
control of the microprocessor 302 drives a graphics display 
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316. A keyboard controller 318 interfaces directly with 
the operator's keyboard 320. An interface port 322 is 
provided for the pointing device 324. 

Operator control of the entire acquisition, analysis 
and display procedure is controlled from the keyboard 320 
and pointing device 324 with feedback, on the graphics 
display 316. One high speed synchronous serial port 326 is 
provided to interface with the 16 channel data acquisition 
system. 1041 Port 326 can be used to send control data to 
the system (e.g., filtering, gain, sampling rate, start/ 
stop acquisition, perform self diagnostics) and to receive 
EC6 data from the system, as well as to download program 
data to the system. Another serial or parallel interface 
port 328 is provided to drive a hard copy output device 330 
for printing desired diagnostic indices. 

Referring now to Fig. 4, the functions performed by 
the system and method of the present invention will now be 
described. As mentioned above, the system and method of 
the present invention quantify linear and nonlinear 
properties of phase and energy components within the 
frequency structure of the ECG from a preselected number of 
leads . Diagnostic indices are then generated from the 
bi spectral data arrays by utilizing predetermined reference 
arrays . The indices are used for the quantification of 
coronary artery disease (CAD), myocardial ischemia, cardiac 
electrical stability, risk of malignant ventricular 
arrhythmia, sit* (si prigin of malignant arrhythmias/ 
extent of malignancy of arrhythmias, degree of * 
antiarrhythmic drug efficacy, neural and humoral inputs to 
the heart, pump function/ ejection fraction, and ongoing # 
organ rejection in cardiac transplant patients. 
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The monitor modulo 402, handles the overall operations 
of the system via integration of data and process 
information from the user interface module 404 , acquisition 
and raw ECG data management module 406, bi spectral 
processing module 408 and diagnostic index derivation 
module 410. A detailed description of the operation of 
modulo 402 will be provided below ,in connection with the 
description of Pig, 5, 

The user interface and display management module 404 
represents the means through which 1 the operator controls 
and interacts with the system during the course of a 
procedure. This includes, but is not limited to, entry of 
information regarding the patient, type of diagnostic 
procedure being carried out; lead and acquisition settings; 
continuous display of acquisition status, lead integrity, 
display of diagnostic indices; and requests for printing 
and, archiving results to disk. Module 404 directly 
interacts with the monitor module 402. The operations 
handled by module 404 can be achieved under one of many 
commercially available environments such as Microsoft's 
Windows. 

The acquisition and raw ECG data management module , 
406, handles all of the raw ECG data checking and 
processing prior to bispoctral analysis. This includes, 
but is not limited to, continuous acquisition of ECG data 
and the verification of the integrity of the data; 
performing QRS detection; performing cross correlation with 
the preselected template to identify suitable complexes; 
extracting suitable complexes from leads of interest in 
preparation for aut obi spectral and crossbi spectral 
processing. Modulo 406 directly interacts with the monitor 
module 402, and a more detailed description of module 406 
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will be provided below in connection with the description 
of Fig. 7. 

The bispectral processing module 408 controls the 
generation of all data arrays measuring dynamic phase and 
energy properties within the ECG. , This information can be 
organized in both aut obi spectral and crossbispectral arrays 
utilizing either an FFT based or parametric based approach. 
The tasks include, but are not limited to, nonlinear/ 
exponential transform of the signal, Fourier 
transformation, the generation of power spectrum, 
autobispectral density, crossbispectral density, 
autobicoherence, crossbi coherence, autobiphase, and 
crossbiphase. Module 408 directly interacts with the 
monitor module 402, and a more detailed description of 
module 408 is provided below in connection with Figs. 9 and 
10. 

The diagnostic index derivation module 410 generates 
the data values utilized in the diagnostic process. The 
task includes, but is not limited to, identifying frequency 
pairs of interest through the use of predetermined clinical 
reference arrays and creating a diagnostic index from the 
values in the bispectral data arrays at the frequency 
locations defined by the reference array. Module 410 
directly interacts with the monitor module 402, and a more 
detailed description of module 410 is provided below in 
connection with Fig. 12; 

^T^t^n^^^bo; Fig . 5 , the operation of th* monitor 
module 402 will now be discussed. In step 502, the data 
arrays used to store the digitized raw ECG and the 
bispectral data of each lead are initialized. The data 
files required for storage and files containing databases 
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required for the computation of diagnostic indices are also 
opened in the initializing step 502. 

( In step 504 the system requests the information 

j required to start the acquisition and diagnostic process 
from the user via the user interface module. This 
information includes , the type of diagnostic procedure to be 
conducted and any operator requested modification to the 
system defaults such as leads used, clinical databases to 
access for diagnostic index computation, filtering, gain 
and sampling rate information for each lead. The type of 
diagnostic procedure as entered by the operator will be 
used by the system to inform the operator of the lead 
placement protocol required, and the type of template to be 
selected (sinus rhythm or ectopic beat) . The system will 
also use the type of diagnostic procedure it will perform 
to select the portion of the complexes to be used as time 
series for bispectral processing as well the type of 
bispectral arrays that need to be computed for use in 
conjunction with the predetermined clinical databases to 
yield a final diagnosis. The process of identifying the 
most effective leads for a particular diagnostic procedure 
and the generation of the clinical databases will be 
discussed later. 

In step 506/ ZCG signals are continuously acquired and 
displayed from the leads needed for the diagnostic 
operation being performed. All channels transmitting 
artif actual data are properly signaled to the operator to 
correct the problem. 'While using the pointing device, the 
operator is requested to select a suitable template against 
which incoming ECG complexes are matched. 

The system, in step, 508, computes the necessary 
autobi spectral and crossbispectral arrays required by the 
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databases for the generation of the diagnostic indices 
requested .by the operator. 

In step 510, the diagnostic indices from all generated 
aut obi spectral and crossbispectral arrays are computed. 
Aut obi spectral density and crossbispectral density clinical 
reference arrays are utilized in these diagnostic index 
computations. In the case of a static one time diagnostic 
determination, such as detection of CAD, the system 
displays/ in step 51 2> the resultant 1 diagnostic indices and 
a projected outcome based on information, from the clinical 
databases. In the case of monitoring an ongoing cardiac 
process such as intraoperative ischemia the index is 
displayed continuously over -time during the course of the 
procedure requiring it . 

In step 514, requested printouts are produced, results 
are stored to disk for archival purposes and all files are 
closed. In step 516, the process is terminated. 

A sample display representation generated by the 
system is shown in Fig. 6. The top section of the screen 
602 is divided into 16 sections 603 each representing the 
region probed by an electrode . Each section will be 
covered by a large w X n 604 if lead fail or artifact was 
detected from the lead corresponding to that section. 

A second portion of the screen 606 can be assigned to 
the continuous display of the diagnostic index if the 
system is being used in monitoring mode. The background of 
that portion is color coded to reflect the possible values 
.allowed for in the range of the selected diagnostic index. 
Tbe most current value of the diagnostic index will dictate 
what color is ^di#plfeyedwi;^fch^backgcou&d^. !! (e.«.9'^=^-- R*d.« 
lowest value to Green « highest value) . This will 
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facilitate the examination of the patient's status at a 
distance. 

A third portion of the scraen 608 can ba assignad to 
displaying one or savaral of tha raw EC6 data laads that is 
baing acquired for processing . This will also provide for 
easy template selection using the pointing device. 

A fourth portion of the screen 610 can be assigned to 
displaying the results of a static diagnostic test if the 
system -is being used in that mode of operation. 

Referring to Pig. 7, the acquisition and raw ECC data 
management module 406 will now be described in greater 
detail. In step 702, the acquisition system 104 is 

programed with requested filtering, gain, sampling rate, 
and lead selection information. 

In step 704, the acquisition system 104 acquires 
continuous ECGdata for all requested leads and transfers 
this data to the host computer 110. The acquisition system 
104 detects lead failures during the acquisition cycle, and 
in step 706, the acquired data is examined for lead failure 
signals and for the presence of artifact. In step 708, 
leads generating fail signals and/ or artif actual data are 
marked for the monitor module 402. 

In step 710, the system requests the operator to 
identify a template from the incoming data stream using the 
pointing device. In step 712, QRS detection is performed 
on the ECG data using any publicly available QRS detection 
program such as the algorithm disclosed by Engeleae at al 
in "A Single Scan Algorithm for QRS -Detection and Feature 
Extraction", IEEE Computers in gajdjfllggy (1979) . Then in 
step 714, the system initiates continuous extraction of 
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complexes that match the template from tha incoming EC6 
data from aach laad using standard cross-correlation 
techniques . 

Depending on the diagnostic test requested by the 
operator and based on information from the clinical 
databases, bispectral processing will be performed on the 
full complex that matches the template or a portion of that 
complex. Also , depending on the test to be performed the 
complex can be a sinus rhythm beat or an ectopic beat, in 
step 716, the auto/cross bispectral time records are; 
continuously generated by extracting a portion suitable for 
the diagnostic test from each matched complex from each of 
the leads required for the test. Bach portion is assigned 
to an Xi(t), where Xi(t) are the individual time series 
records provided for autobi spectral processing. Also in 
step 716, Yi(t) is set to X±+i (t) (the successive complex 
in the same lead) where Yj.(t) are the time series records 
(in addition to Xi (t) ) required for crossbi spectral 
processing within the same lead. Yi(t) can also assume the 
value of a corresponding X±(t) from another lead providing 
for crossbi spectral analysis between two leads. It should 
be noted that for autobi spectral analysis Yi(t) is set to 
equal X±(t) and in all cases the index i denotes the record 
number from 1 to k, and where k is the number of bispectral 
averages computed. 

Once the acquisition process is initiated the program 
* will continuously pass to the monitor module (via interrupt 
handlers) raw time series data until the diagnostic test is 
completed. The program returns control to the monitor 
module 402 in step 718. 

How^,turning^to^Eigi.^ Ai . Fig. a (a) shows a sample PQRST 
template; Fig. 8(b) shows the QRST portion extracted for 
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processing; Fig. 8 (c) shows the terminal part of tha QRS 
extracted for processing; and Fig. 8(d) shows tha ST 
segment with the T wave extracted for processing. The axes 
in Pigs. 8(a) - 8(d) are not calibrated and the figure ara 
for illustration purposes only. 

Referring now to Pig. 9, the frequency domain based 
procedures for producing the autob} spectrum or tha 
crossbispactrum will now be discussed, in step 902, the 
system checks whether the computation to be performed is an 
autobispectral or crossbispectrali computation. 
Autobispectral analysis is a special case of 
croasbispectral analysis and therefore different rules of 
symmetry apply. 

i 

la step 904, the system sets the following aymmetries 
in order to proceed with autobispectral computation: 

f l + f 2 < N/2 

where N is the number of samples in the time series to be 
processed. N will depend on the length of the portion of 
the complex (in msec) and tha sampling rata. Por illustra- 
tion purposes we will utilize the QRS portion of the 
selected complex with length of 160 msec. Since the sam- 
pling rata is 1024 samples/sac N will be equal to 160 
samples or 1024 Hz. 

0 ; < f 2 < f 1 :; 
Xi(t) - Yi(t) — > Xi(f) - Yi(f) 
where fi and f 2 (also referred to aa n and P 2 or Frequency 
1 and Frequency 2) denote the frequency pairs over which 
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bispectral computation will be carried out, Xi {t) and Yi(t) 
denote the individual time series records used for 
bispectral computation, Xi<f) and Y i( f) denote the Fourier 
transform of the time series records and i denotes the 
record number- which! in thia embodiment ranges from 1 to k. 

i * ■ . . ., ' •■ " .■ •' ' ; ' *' . ' ■ ; 

In step 906, the following symmetries are adhered to 
for crossbispectral analysis: 

fi + f2 < N/2 

1 i - 

< 0 < f i < N/2 

0 < f 2 < N/2 
-2f 2 < fi 

Xi(t) * Y A (t) > Xi(f) * Yi(f) 

where all variables represent the same values as they do 
for aut obi spectral analysis, except that for 
crossbispectral analysis Xi(t) and Yi(t> represent in- 
dividually derived time series records where Yi(t) - 
Xi+i<t) from the same lead for same lead crossbispectral 
analysis, or Yi(t) - Xi(t) from a different lead providing 
for inter-lead crossbispectral analysis. 

In step 908, the Discrete Fourier transform (DFT) 
Xi(f) and Yi(f) of each record of the k selected records, 
is computed using a standard IEEE library routine or any 
other publicly available routine in step 908. 

In step 910, the power sp«ctra P x i(f) and F y i (f ) of 
each record of the k selected records is computed by 
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squaring the magnitudes of each .lament of the Pouri.r 
transform Xi<f) and Yi(f) respectively. 

The system computes th. average complex triple product 
xn step 9 12 by utilizing the following eguations where 
bc i (fx / f 2) ia aa individual ceapl-x ^ 

boi(fi,* 2 ) . VxiCfi) * Yi(f 2) * yl(f 1+f2) 
where yl(f 1+ f 2 , ia the complex conjugate of Y i( f 1+ f 2) . 

1 ■ . 

BC(fi,f 2 ) 'm i SbCi(f lr f 2) 

i = 1 

. The average real triple product is computed in step 

indivL'x "V 0lle " ia * br i( *i,f 2 , 
individual real triple product fron, one record and 

BR<fi,f 2 ) is the average real triple product over all 
records: 

bri<fr,f 2 ) - Pxi(fi) * P yi (f 2 , * P yi (f 1+ f 2) 

Jk 

BR(fi,f 2 ) »i £ brKf!/^) 

i = i ' 

In step : 916, the auto/crossbispectral density array 
<BD<fi,f 2 )) is computed using the following equation: 

BD(fi,f 2 ) - |Be<fi,f 2 )| 
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In step 918, the system computes the auto/crossbiphase 
array «p(fi,f2)) using the following equation: 

<P<fi,f2) « tan" 1 [ln(BC(f 1 ,f 2 ))/Re(BC(f 1 ,f2))] 
0 < <p < 2tc (radians) 

In step 920 , the system computes the auto/ 
crossbicoherence array (R (f r , £ 2 ) ) using the following 
equation: i 



R(fl,f2) = BD (fi,f 2 )/[BR(,fi,f 2 H 1/2 
0 < R < 1 

! t 

In step 922, the system returns the requested 
auto/cross bispectral density, average real triple product, 
bicoherence, and biphase arrays to the monitor module 402 . 

Now turning to Fig. 10, a parametric based method for 
producing the autbbi spectrum and the crossbispectrum will 
now be described. In steps 1002, 1004, and 1006 the system 
sets the symmetries and time series records in the same 
manner as described above in steps 902, 904, and 906, 
respectively. The power spectra of x±(t) and Yj.(t) are 
estimated in steps 1008, 1010, and 1012. This estimation 
method includes two major stages, the autoregressive (AR) 
model order selection and the power spectrum computation 
for Xi (t) and Yi(t) . In step 1008/ the system computes two 
sequences of autocorrelations; {R 2x (m) } and {R2y(m) } using 
the following equation. 

h w-imi 

R2M " M~7tT S X ^i(t)zi(t+m), m = 0, 1, L, z^X±,Yl 



WO 91/02484 



PCT/US90/01385 



-27- 

where M is the number of records (Jc in our case) , and N is 
the number of samples per record (160 in our case) . L is 
much greater than the possible AR filter order (we choose 
50). 

I. . 

The Final Prediction Errors, FPEx(m) and FPEy (m) are 
calculated for all orders, m - 1, 2, ... t; by performing 
a Levinson recursion function on each autocorrelation 
sequence in step 1010 in order to find the order of the AR 
filter. The locations, Qx and Qy, of the minimum of 
FPEx(m) and FPEy(m) respectively are chosen to be the 
orders of the AR filters of power spectra of Xi (t) and 
Yi(t) respectively, i.e., 

FPEx(Qx) = min {FPExCm)} and FPEy{Qr) = min {FPEy(m)} 

Once the orders of the AR filters, for power Spectra 
are chosen, the autocorrelation sequences, (R2x(m) } and 
(R2Y(m)), are entered into Levinson recursion with order Qx 
and Qy, respectively, instead of L. The coefficients, 
(cix, i«6, 1, ...,Qx) and {ciy, i » 0,1, ... ,Qy>, obtained 
from the recursion are the coefficients of the AR filters 
for power spectra of Xi(t) and Yi(t) respectively. Then, 
in step 1012, the power spectrum, P*(f), is computed as the 

prediction error (of) divided by square of the magnitude of 
the Fourier transform of the coefficients, i.e., 

■» ■ J. 

Pzif) = j£ 2 -, z = X,Y. 

I 1 + 2 ciz e-X* fi I 2 

i -. 1 
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The system estimates the auto/cross bispectrum in 
steps 1014, 1016, and 1018. The estimation process 
includes two major stages: the order selection and 
bispactrum computation. In step 1014, two sequences of 
third-order moments, {R3x(T)} and {R 3 y(t)> are computed 
using: the following equation. 



*3z(x) = — - £ J *i(t)fi(t+c), z = x, Y, and T = —L, 



i=l t=sx 

where si -max (1,1-t), a? » min (N, N-x), and L is much 
greater than the possible AR filter orders (e.g. SO). 

In step 1016, two super matrices T x and Ty are formed 
as follows. 



2V = 



R3z(rL-l) R3z(-L) j J? 3z (-l) 

R3z(~2L) ,R 3z (-2L+l) ... R3 z (rL) 



X, Y 



From the assumption we made about the AR filter of 
bispactrum, the orders Ox and Oy of the AR filters of 
bispectra of X±(t) and Y± <t) are the ranks of the super 
natricee Tx and Tr. Therefore, Ox and Che *re chosen by 
using singular value decomposition . Having found the 
orders, we obtain the coefficients of the AR filters of 
bispectra by solving the following linear system of 
aquations: 



f*3*(0) R3x(1) 
*3*<-l) R3z(0) 



R3z(-0z) R2z(r0z+D 

V • ■■■ - ■ 



RZziOz) A 
R3z{0z-1) 



*3z(0) 



A 



bl; 



i>o zZ 



■ / 



0 



0 



X, Y. 
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vhere th. skewness <p z) and ^ coefficients (b, 2 

bo z z) , x « X, Y, can be obtained by solving the linear ays- 

tern of equations . 

The auto/cross bispectrum of x i( t) and y i( t) are 
computed in step 1018 as the cubic root of the tripi. 
product of the skewnesses divided by the triple 

product of the Fourier transforms of the AR filter 
coefficients (H z (f)), i.«./ i 



BC{fi,f 2 ) « (My0y) V V^ 1 )HKf2)Hkfl+f2) 
H z if .) « 1 + £ bi- , r « x,y. 



and BR (fx, fj) is the real triple product for that same lead: 
BR(fi,f2) - Px(fl) * Py(f 2 ) * P y (f 1+ £ 2 ) 



After obtaining power spectrum and auto/cross 
bispectrum, the system computes the bispectral density 
array, the biphase, and the bicoherence in step 1020 the 
same way as in steps 916, 918, 920, In step 1022, the 
aystem returns to the monitor module 402 the requested 
bispectral density,, biphase, and bicoherence arrays. 

For illustration purposes Pig. 11 contains sample 
autobispectral arrays of the W comply showing frequency 
pairs 0 < fi < 512 Hz, and 0 < f 2 < 256 Hx. A bispectral 
density array, is shown in Pig, 11(a) where the z axis 
represents the magnitude in decibel. (db) r of Uncoupling 
interaction between all appropriate frequency pairs fl and 
f 2 . Recall that the frequency pairing scheme must adhere 
to symmetry rule: 
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f 1 + f 2 < N/2 

where H » 1024 Hz in this cast. A biphase array is shown 
in Pig. 11 (b) whore tha Z axis represents the phase in 
radians of the coupling interaction between all appropriate 
frequency pairs fi and f 2 . A bi coherence array is shown in 
Pig. 11(c) where the Z axis represents the normalized mag- 
nitude in percent (%) of the coupling interaction between 
all appropriate frequency pairs fi and £2. 

•1 1 

Referrihg now to Fig. 12 a mora detailed outline of 
tha diagnostic index generation module 410 is shown, in 
step 1202, the software identifies the type of diagnostic 
test in progress. The possible options include but are not 
limited to: 

1. detection and quantification of coronary artery 
disease 

2. detection and monitoring of myocardial ischemia 
3 • dat0ctiof ? of re P* rf us^ion status post throm- 
bolytic therapy 

4. assessment of coronary artery restenosis after 
successful percutaneous transluminal coronary an- 
gioplasty 

5. one- time or continuous monitoring of chances in 
cardiac electrical stability whether due to drugs, 
heart disease, acute infarction or neural 'factors. v '" 

6. one-time or continuous assessment of antiarrhyth- 

mic dxrujj ef fe^ts ,,c^c^ 
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7. assessment of the extent of malignancy of car* 
diac arrhythmias. 

8 • differentiation of sustained ventricular 
tachycardia from wide complex supraventricular 
tachycardia. ' 

9. one time or continuous monitoring of neural in- 
puts to the heart, 

i 

\0. evaluation of pump function and ejection frac- 
tion 

11. identifying ongoing organ rejection in cardiac 
transplant patiepts - 

In step 1204 , the system retrieves the type of 
bispectral array to use in the diagnostic index computation 
as well the lead(s) of origin. We have 16 possible EC6 
leads and each lead could have 4 aut obi spectral arrays and 
4 crossbispectral arrays for any portion of the complexes 
matching the template selected. These arrays do not 
include inter-lead crossbispectral arrays. The combination 
possibilities can lead to a very large number of 
computations associated with each diagnostic index as well 
as a substantial reference clinical database. Statistical 
methods to identify leads and bispectral arrays with the 
greatest diagnostic fidelity for each test are used to make 
this system practical. Such statistical methods will be 
discussed in greater detail below. 

In stop 1206, the appropriate reference array is 
retrieved from resident memory or from disk. Each 
reference array will contain the locations of the frequency 
pairs which are most sensitive to the diagnostic test in 
progress. In step 1208, all data points in the bispectral 
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array at the locations identified by tha retrieved 
reference array ara addad togathar for a singla value 

index. A counter (NP) of tha total number of points addad 
is leapt. In step 1210, tha singla valua indax is divided 
by NP to obtain the diagnostic index. In step 1212, the 
program returns to the monitor module 402 . . 

The predetermined clinical reference arrays referred 
to above are critical to the device's ability to achieve 
clinically relevant diagnostic efficacy, and the process 
adopted for generating these clinical reference arrays will 
now be described. Since the total number of possible 
diagnostic applications will require many reference arrays, 
only two types of statistical approaches will be discussed 
in detail. All other reference arrays are acquired in a 
similar fashion utilizing a wide number of clinically 
appropriate statistical approaches. For illustration 
purposes the generation of tha bicoharence reference array 
for detection of coronary artery disease with orthogonal 
lead X, assessment of risk to VT/VF with orthogonal lead Z, 
and assessment of ejection fraction with orthogonal lead X 
will be discussed in the following section. 

In a first study raw ECG signals were acquired from 
two groups of subjects: 

a) young, medically healthy normal volunteers (n«28) 

b) older patients with normal conventional ECG's 
and angiography-proven coronary atherosclerosis 
<n-16) 

The acquisition procedure descrJU^d previously is 
followed: 

- Band pass filter 0.05-512 Ha 

— Sampling rate- 1024- sample/sec 
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- Number of bispectral avrag., performed k-100 
For all subjects thr« minut. . of artifact free ECG 
data ware acquis. An autobicobaranca array i 3 generated 
from 100 QRS complexes from ortbogonal lead x for £ 
.ubject.. Tb. arrays aragroupad in 2 sat. of arrays, the 
first representing the young normals and tba second 
raprasanting tba coronary artary disease subjects. 

t 

A pairad Student's t test is Mrf Am ^ - . 

*■ *« parxormod on aacb of 1640 

data points, comparing tba first and sacond sats of arrava 
The resulting 1640 t values ara adored in a two dimensional 
array identical in structure to tbat of tba bicobarence 
array. 

Pach t value from t array (T(fl,f2>> is tested for 
significance based on tba number of degrees of freedom 
Wbere tbe degrees of freedom are equal to tba total number 
of subjects-!. All t value, not meeting tb. required 
significance level are set to 0. m tbe preferred 
ambodiment all locations vitb a t value not corresponding 
to a p < 0.05 are set to 0. 

Tbe application of tba above conditions has tbe effect 
of identifying all of tbe frequency pair locations tbat are 
significantly different in coronary artery patient, when 
compared to young normal volunteer.. Suppression of all 
other frequency pair, allow, easier inspection of tbe most 
sensitive regions. 

Turning now to Pig. 13, Pig. 13(a) .bow. tb. mean 
autobicobaranca array for ortbogonal lead X for tbe normal 
subject.. Pig. 13(b) show, the mean autobi coherence array 
for orthogonal lead Y f or the CAT) subject.. Pig. 13(e) 
show, the t array with all t value, not meeting p < o.OS 
aet to 0. 
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The -next stop involves sorting the t array for the 
most sensitive ensemble of frequency pair locations. In a 
preferred embodiment this would consist of the top. 25% of 
all significant t values. The locations -fl/£2 of the most 

5 significant t values will be used to generate a diagnostic 
index for each subject in the process described above. 
Fig* 14 (a) shows sample .indices for normal subjects 
and Fig,* 3.4 (b) ' shows sample indices for CAD 
subjects. The coded filename of the subject tested 

10 precedes each index. 

i , 

The final step is to identify a cutoff value for the 
diagnostic index above which subjects are normal and all 
values below will indicate the presence of CAD . This 
cutoff should be optimized to yield the best sensitivity 
15 and specificity. 

♦ 

. ; tr ue positive tests 

Sensitivity = all tzue positives 

Specificity = true ne< f ative tests 
speciticxzy all true negatives 



In this particular case if wtj chost a entoff of 20.0 
the sensitivity and specificity will ba 100% (16/16) and 
20 96.4% (27/28) respectively. 

The .above,, statistical steps are repeated for all the 
possible bi spectral arrays that can ba generated (as 
dascribad above) from all ECG leads acquired. To rank 
ordar tha arrays with respect to diagnostic efficacy * 
25 second prospective study is conducted. The conditions ? 
under which the study, is conducted are identical to those 
of the first except that: a) the frequency pair locations 
of interest have already been identified and are now 
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followed prospectively and b) the size of the study group 
is now sufficiently large so that sample variation of 
bxspectral arrays more closely approximates the true 
variance within the population undergoing the intervention 
or suffering from the disease. 

i - 

Thus for the example of coronary artery diagnosis, 
patients undergoing coronary angiography would each have, a 
recording taken. The diagnostic index for each type of 
bispectral array from each lead ,©r combination of leads 
would be calculated. These indices would then be compared 
to the results of angiography. The sensitivities and 
specificities for the detection of coronary disease would 
then be calculated for each bispectral array. The best 
array and lead system is then chosen as the default setting 
for the diagnostic procedure under consideration. The 
frequency Ideations identified by t values of this array 
will be used by the diagnostic index derivation module 410. 

The above statistical approach may also be used to 
identify patients who are at risk of malignant arrhythmias . 
in a second study two new groups of patients are used: 

a) The control group consists of patients with pre- 
vious anterior or inferior MX without a history, of 
malignant ventricular arrhythmias (n-37) . 

b) The arrhythmia group consists of patients with a 
history of previous anterior or inferior MX and a 
history of sustained ventricular tachycardia or out- 
of -hospital cardiac arrest (VT/VP) (n*12) . 

An autobi coherence array is generated from 100 QRS 
complexes from orthogonal lead Z for all subjects . The 
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arrays ara grouped in 2 seta of array a, the first 
representing tha MX subjects and the aecond representing 
tha VT/VT subjects. Aftar conducting tha t test described 
aarliar tha statistical data shown in rig. 15 was 
5 ganaratad. Fig. 15(a) shows tha mean autobicoharanca array 
for orthogonal lead Z for tha Ml subjacts. Pig. 15(b) 
shows tha mean autobicoharanca array . for orthogonal lead Z 
for tha VT/VT subjacts. Pig. 15(c) shows tha t array with 
all t values not meeting p < 0.05 sat to 0. 

10 tha t array ganaratad is than usad to product 

diagnostic indices for both groups as Shown aarliar for tha 
CAD study. Tha whola statistical process is followed 
through including prospective studies to idantify tha bast 
bispactral array and its corresponding laad for tha 
15 idantification of patients who are at risk of malignant 
arrhythmias. Fig. 16(a) shows sample indices for 
MI subjects and Fig. 16(b) shows sample indices for 
the VT/VF subjects. The coded filename of the subject 
tasted pracadas each index. In thia particular case if we 
20 chose a cutoff of 26.0, the sensitivity and specificity 
will be 91. 6% (11/12) and 86.4% (32/37) respectively. 

This invention is not limited to the use of the t test 
and many other statistical ranking tests might be used when 
appropriate. For example if the diagnostic procedure 
requires the measurement of a continuous variable such as 
ejection fraction (EF) regression analysis is mora suitably 
used, ^".j . h 

The following example will describe the use of the 
present invention to measure a continuous variable. In 
this example, raw ECG signals were acquired from a large 
number of patients (n*247) with known EF through 
ventriculography, the autobicoharanca array from 100 QRS 

SUBSTITUTE SHEET 
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co» P l«„ for l..d x »,. eonpufd. Eaoh on. or th. 1640 
bicoh.r en =a point, w„ corralafd with it, oo«.,pondin9 EF 
vain, aoro., .11 p.ti.nt.. Th. o«too«. i. a t«, 
di».n,ion,l ,rr.y ,i^ llr ia 3t ructur, to that o* th. 
bieoherene. ,rray and wh.r. th. vain. , t Mch .„ „. 
r. P r.,.nt. th. r r.,r.„io„ oo.«ioi.»t b.t«..n th. 
bxooh.r.nc. at thi, Ir.qn.ncy pair Md . imatttm tractton 
Th. r oo.«i=i.„ t » ilr rang . fro- ■ M 

•hown inpi g . 17 . Pig . 17(m) ,ho». th. po.iti«7^:.: 
b.t„„ „ „ d Mcoh .„ nc . , t th . co „ Mpo £r> ^ 
pa«.,-,nd rig. 17(b) ,ho». th. n.g.tiv r «i„... °* 

In addition to quantifying the presence and axtant of 
coronary artery disease (CAD) , myocardial ischa^^ cardiac 
•metrical stability, risfc of -aligns vantric^ar ^ 

ext!»T?' W 0t ° rigin ° £ arrhythmias, 

axtant of malignancy of arrhythmiaa, degree of 

antiarrhythmic drug efficacy, neural and humoral inputs to 
tha heart, pump function/ auction fraction, and ongoino 

IT* TT l9n ^ C3rdiaC tranS * lant Pliant,, the sys L 
and mathod of tha prasent invention may also ba u.ad to 
assess a myriad of cardiac phenomena baead oh the 
acquisition and processing of 2C6 signals into variou. 
bispectral arrays which ara then compared to appropriata 
reference arrays. , 

Although power spectral analytic techniques in the 
frequency domain hav. been applied to tha ECG signal as 

^n"^ 10 ^ BaCkgrOUn d ^ orde^ spectra 

whether by pft or by parametric approaches hav. never been 

so applied. No bispactral technique has aver bean 

demonstrated to be useful for any diagnostic purpbsa. 

Other techniques for the quantification of coronary artery 

disease, PTCA restenosis or the detection of cardiac 

•lactrical instability, aapecially aa it changes with 
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antiarrhythmic drug administration, remain qualitative , 
static and limited in their overall utility and acceptance 
in practice. Specifically, the system and method of the 
present invention uses various bispectral values to measure 
dynamic frequency structure [higher order phase-locking] 
across all frequency pairs in a frequency range Ignored by 
those knowledgeable in the art and Uses various alterations 
in these bispectral parameters at a limited number of 
frequency locations as an index of physiological 
perturbation . The system and method utilizes various 
bispectral arrays of transformed ZCG signals of defined 
clinical populations to define the locations of the subset 
of frequencies used to calculate this index. Reference 
clinical arrays are further utilized to assess the meaning 

of this index and to measure the significance of deviations 

-■ • ■■ ■ it . 

of this index from normality. This allows the quantitative 
gauging of the disturbances in cardiac function, whether 
due to coronary disease, electrical instability, restenosis 
after PTCA, drugs or ischemia for any particular ECG lead 
position. The invention disclosed here also defines the 
graphic display of the diagnostic index, whether on video 
screen or on paper, whether in real-time or in digital 
archive... , 

While the foregoing invention hat btan described with 
reference to its preferred embodiments, various alterations 
and modifications will occur to those skilled in the art. 
For example, the ECG may carry diagnostic information at 
frequencies much higher then the cutoff frequency of 512 
Hz; ^she-uae^ra^^uch^h^gb-f requene^.4c««r.energy,d6B^ohehts : - 
of the ECG waveform by the system and method described 
above is intended to fall within the scope of the current 
invention. All auch alterations and modifications are 
intended to fall within the scope of the appended claims. 
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What is claimed is: 

1. A method of noninvasively detecting cardiac 
phenomena comprising the steps of: 

acquiring electrocardiographic signals from a 
body surface of a subject being analyzed through a surface 
electrode; 

i 

determining a waveform template representing a 
QRST waveform in said subject; 

selecting a number of signals that match said 
waveform template and storing said signals for processing; 

generating bispectral values from selected 
signals that can be used to analyze cardiac phenomena. 

2. The method of noninvasively detecting cardiac 
phenomena of claim 1 wherein said generated bispectral 
values are autobi spectral density values. 

3. The method of noninvasively detecting cardiac 
phenomena of claim 1 wherein said generated bispectral 

values are autobicoherenco values. 

4. The method of noninvasively detecting cardiac 
phenomena of claim 1 wherein said generated bispectral 
values are autobiphase values. 

5. The method of noninvasively detecting cardiac 
phenomena of claim 1 wherein said step of selecting a 
number of signals that match said waveform template 
comprises the step of selecting pairs of successive time 
series records that are used in the generation of 
bispectral values. 
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6. The method of noninvasively detecting cardiac 
phenomena of claim 5 wherein said generated bispectral 
values artf crossbispectral density values. 

7. The method of noninvasively detecting cardiac 
phenomena of claims 5 wherein said generated bispectral 

values are crossbicoherence values. 

1 i . ...... 

8. The method of noninvasively detecting cardiac 
phenomena of claim 5 wherein said generated bispectral 
values are crossbiphase values. 

9 « The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the step of using 
said bispectral values to detect a presence and extent of 
coronary artery disease in said subject. 

10. The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the step of using 
said bispectral values to detect a presentee and extent of 
myocardial ischemia in said subject. 

11 . The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the step of using 
said bispectral values to analyze cardiac electrical 
stability in said subject. 

12. The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the step of using 
said bispectral values to analyze a risk of malignant 
ventricular arrhythmia in said subject. 
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13. The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the step of using 
said bispectral values to locate a site of origin of 

1 malignant arrhythmias in said subject. 

14. The method of noninvasively detecting cardiac 
.- •• i 

phenomena of claim 1 further comprising the step of using 

said bispectral values to analyze an extent of malignancy 
of arrhythmias in said subject. 

i ■* ■ .* . • 

15. ^The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the step of using 
said bispectral values to evaluate a degree of 
antiarrhythmic drug efficacy in said subject. 

16. The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the step of using 
said bispectral values to examine neural and humoral inputs 
to the heart of said subject. 

17. The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the step of using 
said bispectral values to evaluate pump function or 
ejection fraction in the heart of said subject. 

18. The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the step of using 
said bispectral values to monitor ongoing organ xe jection 
wherein said subject is a cardiac transplant patient. 

19. The method of noninvasively identifying and 
distinguishing wide-complex supraventricular tachycardia 
from sustained ventricular tachycardia. 
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20. The method of noninvasively detecting cardiac 
phenomena of claim 1 wherein said step of generating 
bispectral values comprises: 

assigning each selected signal matching said 
waveform template to a time series record; 

computing an average complex triple product using 
Fourier transforms of said time series records; 

computing biphase values as, the inverse tangent 
of the quotient of an imaginary part of the average complex 
triple product divided by a real part of the average 
complex triple product . .-■»«.. 

21. The method of noninvasively detecting cardiac 
phenomena of claim 20 wherein said biphase values are 
crossbiphase values. 

22 . The method of noninvasively detecting cardiac 
phenomena of claim 20 wherein said biphase values are 
autobiphase values. 
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23. The method of noninvasively detecting cardiac 
phenomena of claim 1 wherein said step of generating 
bispectral values comprises: 

v assigning each selected signal matching said 
waveform template to a time series record; 

determining Fourier transforms of said time, 
series records; 

estimating bi spectrum of said time records by 
computing sequences of third-order moments, forming super 
matrices using said third-order moments, determining orders 
of autoregressive model filters of the bi spectrum of the 
Fourier transforms of said time records from ranks of said 
super matrices, obtaining the skevness and autoregressive 
filter coefficients by solving a set of equations based on 
said third order moments and said' orders, and computing 
said bispectrum from said skewness and autoregressive 
filter coefficients. 

« 

24. The method of noninvasively detecting cardiac 
phenomena of claim 23 further comprising the step of 
estimating power spectra of said Fourier transforms, said 
step of estimating power spectra comprising the steps of: 

computing two sequences of autocorrelations and 
performing a recursion function on each autocorrelation 
sequence in order to determine final prediction errors, 
minimums of said final prediction mxrora being orders of 
the autoregression filters of said power spectra; 

performing a recursion function on said 
autocorrelation sequences using said orders of the 
autoregression filters in order to generate coefficients of 
said autoregression filters; 

computing the power spectra as the final 
prediction error divided by a square of the magnitude of 
the Fourier transform of said coefficients. 
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25. The method of noninvasively detecting cardiac 
phenomena of claim 23 wherein said step of computing 
bispectrum of said Fourier transforms of said time series 

'' records comprises dividing the cubic root of the product of 
* three skevnesses by a triple product of the Fourier 
transforms of the AR filter coefficients- * 

26. The method of noninvasively detecting cardiac 
phenomena of claim 23 further comprising the step of 
computing an array of bispectral densities by determining 
the absolute value for each bispectrum value. 

27. The method of noninvasively detecting cardiac 
phenomena of claim 23 further comprising the step of 
computing an array of biphase values as the inverse tangent 
of the quotient of an imaginary part of the bispectrum 
divided by a real part of the bispectrum. 

28. The method of noninvasively detecting cardiac 
phenomena of claim 24 further comprising the steps of: 

computing a real triple product as the product of 
three power spectra values; 

computing an array of bicoharence values as the 
quotient of bispectral density divided by a square root of 
said real triple product . 

29. The method of ndhinvaSively detecting cardiac 
phenomena of claim 1 wherein said bispectral values are 
crossbi spectral values generated from consecutive signals 
from one surface electrode. 
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30. The method of noninvasively detecting cardiac 
phenomena of claim 1 wherein said bispectral values are 
crossbispectral values generated from signals obtained from 
two surface electrodes. 

31. The method of noninvasively detecting cardiac 
phenomena of claim 1 wherein said template is a sinus 
rhythm template. 

32. The method of noninvasively detecting cardiac 
phenomena of claim 1 wherein said template is an ectopic 
beat template. 

33. The method of noninvasively detecting cardiac 
phenomena of claim 1 further comprising the steps of 
generating clinical reference arrays for use in identifying 
particular bispectral values that are sensitive to the 
cardiac phenomena being analyzed. 

34. The method of noninvasively detecting cardiac 
phenomena of claim 33 furthir comprising the step of rank 
ordering said clinical reference arrays with respect to the 
diagnostic efficacy of the arrays for the cardiac phenomena 
being analyzed. 

■ 

35. The method of noninvasively detecting cardiac 
phenomena of claim 33 wherein values in said clinical 
reference arrays that are not significant are suppressed to 
facilitate inspection of the most sensitive values . 
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36. The method of noninvasively detecting cardiac 
phenomena of claim 33 further comprising the steps of 
determining a diagnostic index from said clinical reference 

f array, generated values above said index indicating a 
normal condition and below said index indicating an 
abnormal condition. 

37. The' method of noninvasively detecting cardiac 
phenomena of claim 36 wherein said diagnostic index is 
optimized to yield the best sensitivity and specificity 
with sensitivity being equal to true positives divided by 
positives and specificity being equal to true negatives 
divided by negatives • 

38 . A system for noninvasively detecting cardiac 
phenomena comprising: 

means for acquiring electrocardiographics signals 
from a surface of a body of a subject being analyzed; 

means for determining a waveform template 
representing a normal QRST waveform in said subject; 

means for selecting a number of signals that 
match said waveform template; 

means for storing said signal portions for 

processing; 

means for generating bi spectral values from 
selected signals that can be used to analyze cardiac 
phenomena* 

39 . The system for noninvasively detecting cardiac 
phenomena of claim 38 wherein said bispectral values are 
crossbiphase values. 
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40. The system for noninvasively detecting cardiac 
phenomena of claim 38 wherein said bi spectral values are 
autobiphase value. 

41. The system for noninvasively detecting cardiac 
phenomena of claim 38 wherein said generated bispectral 
values are autobispectral density values. 

42. The system for noninvasively detecting cardiac 
phenomena of claim 38 wherein said generated bispectral 
values are crossbispectral density values. 

43. The system for noninvasively detecting cardiac 
phenomena of claim 38 wherein said generated bispectral 
values are aut obi coherence values . 

44. The system for noninvasively detecting cardiac 
phenomena of claim 38 wherein said generated bispectral 
values are crossbi coherence values. 

45. The system for noninvasively detecting cardiac 
phenomena of claim 1 wherein said means for acquiring 
electrocardiographic signals is a plurality of surface 
electrodes and wherein said bispectral values are 
crossbispectral values generated from consecutive signals 
from one surface electrode. 

46. The system for noninvasively detecting cardiac 
phenomena of claim 38, wherein said means for acquiring 
electrocardiographic signals is a plurality of surface 
electrodes and wherein said bispectral values are 
crossbispectral values generated from signals obtained from 
two of said surface electrodes. 
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47 . The system for noninvasively detecting cardiac 
phenomena of claim 38 wherein said template is an ectopic 
baat tamplata. 

48. The system for noninvasively detecting cardiac 
phenomena of claim 38 further comprising means for 
generating clinical reference arrays for use in identifying 
particular bispectral values that are sensitive to the 
cardiac phenomena being analyzed. 

49. The system for noninvasively detecting cardiac 
phenomena of claim 38 further comprising means for 
automatically rank ordering said clinical reference arrays 
with respect to diagnostic efficacy of the arrays for the 
cardiac phenomena being analyzed. 

50. The system for noninvasively detecting cardiac 
phenomena of claim 38 further comprising means for 
suppressing values in said clinical reference arrays that 

are not significant in order to facilitate inspection of 
the most sensitive values. 

51. the system for noninvasively detecting cardiac 
phenomena of claim 38 further comprising means for 
determining a diagnostic index from said clinical reference 
arrays, generated values above said index indicating a 
normal condition and values below said index indicating an 
abnormal condition. 
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S2. The system for noninvasively detecting cardiac 
phenomena of claim 38 wherein said diagnostic indaac is 
optimized to yield the best sensitivity and specificity 
with sensitivity being equal to true positive tests divided 
by all true positives and specificity being equal to true 
negative tests divided by all true negatives. 

• V ■ 
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AMENDED CLAIMS 
[received by the International Bureau 
on 21 December 1990 (21.12.90); 
original claims 1,9-19,38,40,45,48-52 amended; 
new claims 53-56 added; other claims unchanged (8 pages)] 

4 

H 1 . A method of noninvasively detecting cardiac 

phenomena comprising the steps of: ' ~ 

acquiring electrocardiographic signals from a 
body surface of a subject being analyzed through a 
surface electrode; 

determining a waveform template representing a 
QRST waveform in said subject; 

selecting a number of signals that match said 
waveform template and storing said signals for 
processing; i , 

generating bispectral values from selected 
signals, and comparing saxd bispectral values to 
reference bispectral values to analyze cardiac 
phenomena. \ 

2. The method of noninvasively detecting 
cardiac phenomena of claim 1 wherein said generated 
bispectral values are autobispectral density values* 

3. The method of noninvasively detecting 
cardiac phenomena of claim 1 wherein said generated 
bispectral values are autobicoherence values* 

4. The method of noninvasively detecting 

cardiac phenomena of claim 1 wherein said generated ^ 
bispectral values are autobiphase values. 

\^ 

5. The method of noninvasively detecting * 
cardiac phenomena of claim 1 wherein said step of 

selecting a number of signals that match said waveform 
template comprises the step of selecting pairs of 
successive time series records that are used in the 
generation of bispectral values. 
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6. The method of noninvasively detecting 
cardiac phenomena of claim 5 wherein said generated 
bispectral values are crossbispectral density values. 

i 

7. The method of noninvasively detecting 
cardiac phenomena of claims 5 wherein said generated 
bispectral values' are crossbicoherence values. 

8. ' The method of noninvasively detecting 
cardiac phenomena of claim 5 wherein said generated 
bispectral values are crossbiphase values. 

9. The method of noninvasively detecting cardiac 
phenomena 1 of claim 1 further comprising the step of 
comparing said bispeqtra! values to reference 
bispectral values to detect a presence and extent of 
coronary artery disease in said subject. 

10. The method of noninvasively detecting 
cardiac phenomena of claim 1 further comprising the 

step of comparing said bispectral values to reference 
bispectral values to detect a presence and extent of 
myocardial ischemia in said subject. 

11. The method of noninvasively detecting 
cardiac phenomena of claim 1 further comprising the 
step of comparing said bispectral values to reference 
bispectral values to analyze cardiac electrical 
stability in said subject. 

12. The method of noninvasively detecting 
cardiac phenomena of claim 1 further comprising the 
step of comparing said bispectral values to reference 
bispectral values to analyze a risk of malignant 
ventricular arrhythmia in said subject. 
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13. The method of noninvasively detecting 
cardiac phenomena of claim 1 further comprising the 
step of comparing said bispectral values to reference 
bispectral values to locate a site of origin of 
malignant arrhythmias in said subject. 

14. The method of noninvasively detecting 
cardiac phenomena of claim 1 further comprising the 
step of comparing said bispectral values to reference 
bispectral values to analyze an extent of malignancy 
of arrhythmias in said subject. 

15. The method of noninvasively detecting 
cardiac phenomena of claim 1 further comprising the 
step of comparing said bispectral values to reference 
bispectral values to evaluate a degree of 
antiarrhythmic drug efficacy in said subject. 

16. The method of noninvasively detecting 
cardiac phenomena of claim 1 further comprising the 
step of comparing said bispectral values to reference 
bispectral values to examine neural and humoral inputs 
to the heart of said subject. 

17. The method of noninvasively detecting 
cardiac phenomena of claim 1 further comprising the 
step of comparing said bispectral values to reference 
bispectral values to evaluate pump function or 

e jection ; f ractipn i^fehe he^rt^ pf said subject. 
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18. The method of noninvasively detecting 
cardiac phenomena of claim I further comprising the 
step of comparing said bispectral values to reference 
bispectral values to monitor ongoing organ rejection 
wherein said subject is a cardiac transplant patient. 

19. The method of noninvasively detecting 
cardiac phenomena of claim 1 further comprising the 
step of identifying and distinguishing wide-complex 
supraventricular tachycardia from sustained 
ventricular tachycardia. 

20. , The method of noninvasively detecting 
cardiac phenomena of claim 1 wherein said step of 
generating bispectral values comprises: 

assigning each selected signal matching said 
waveform template to a time series record; 

computing an average complex triple product 
using Fourier transforms of said time series records? 

computing biphase values as the inverse 
tangent of the quotient of an imaginary part of the 
average complex triple product divided by a real patt 
of the average complex triple product. 

21. The method of noninvasively detecting 
cardiac phenomena of claim 20 wherein said biphase 
values are crossbiphase values . 

22. The method of noninvasively detecting 
cardiac phenomena of claim 20 wherein said biphase 
values are autobiphase values. 
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36. The method of noninvasively detecting 
cardiac phenomena of claim 33 further comprising the 
steps of determining a diagnostic index from said 
clinical reference array, generated values above said 
index indicating a normal condition and below said - 
index indicating an abnormal condition* 

37. The method of noninvasively detecting 
cardiac phenomena of claim 36 wherein said* diagnostic 
index^is optimized to yield the best sensitivity and 
specificity with sensitivity being equal to true 
positives divided by positives and specificity being 
equal to true negatives divided by negatives. 

38. A system for noninvasively detecting cardiac 
phenomena comprising,: 

means for acquiring electrocardiographic signals 
from a surface of a body of a subject being analyzed; 

means for determining a waveform template 
representing a normal QRST waveform in said subject; 

means for selecting a number of signals that 
match said waveform template; 

means for storing said selected signals for 
processing; 

means for generating bispectral values from said 
selected signals and for comparing said bispectral 
values to reference bispectral values to detect and 
analyze cardiac phenomena. 

39 4 . The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said bispectral 
values are crossbiphase values. 
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40. The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said bispectral 
values are autobiphase values, 

41 • The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said generated 
bispectral values are autobispectral density values. 

42. The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said generated 
bispectral values are crossbispectral density values. 

43. The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said generated 
bispectral values are autobicoherence values. 

44. The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said generated 
bispectral values are crossbicoherence values. 

45. The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said means for 
acquiring electrocardiographic signals is a plurality 
of surface electrodes and wherein said bispectral 
values are crossbispectral values generated from 
consecutive signals from at least one of said 
plurality of said surface electrodes. 

46. The, system for noninvasively detecting 
cardiac phenomena of claim 38 , wherein said means for 
acquiring electrocardiographic signals is a plurality 
of surface electrodes and wherein said bispectral 
values are crossbispectral values generated from 
signals obtained from two of said surface electrodes. 
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47. The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said template is 
an ectopic beat template. 

48. The system for noninvasively detecting 
cardiac phenomena of claim 38 further comprising means 
for generating clinical reference arrays for use in 
identifying particular bispectral values that are 
sensitive to the cardiac phenomena being detected and 
analyzed. 

49. The system for noninvasively detecting 
cardiac phenomena of claim i38 further comprising means 
for automatically rank ordering clinical reference 
arrays of said reference -bispectral values with 
reS p ec t to diagnostic efficacy of the arrays for the 
cardiac phenomena being detected and analyzed. 

50. The system for noninvasively detecting 
cardiac phenomena of claim 38 further comprising means 
for suppressing values in clinical reference arrays of 
said reference bispectral values that are not 
significant in order to facilitate inspection of 
values most likely to indicate the presence of the 
cardiac phenomena. 

51. The system for noninvasively detecting 
cardiac phenomena of claim 38 further comprising means 
for determining a diagnostic index from clinical 
reference arrays of said reference bispectral values r 
generated values above said index indicating a normal 
condition and values below said index indicating an 
abnormal condition. 
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52. The system for noninvasively detecting 
cardiac phenomena of claim 51 wherein said diagnostic 
index is optimized to yield the best sensitivity and 
specificity with sensitivity being equal to true 
positive tests divided by all true positives and 
specificity being equal to true negative tests divided 
by all true negatives. 

■ i . 

1 ' ■ ■ 

53. The method of noninvasively detecting 
cardiao phenomena of claim 1 Wherein said bispectral 
values are generated by computing the Fourier 
transform of the third order autocorrelation function 
of said filtered signals. 

54. The method of noninvasively detecting 
cardiac phenomena of claim 1 wherein said bispectral 
values are' generated by computing the Fourier 
transform of the third order crosscorrelation function 
of said filtered signals. 

55. The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said means for 
generating bispectral values comprises a means for 
computing the Fourier transform of the third order 
autocorrelation function of said filtered signals. 

56. The system for noninvasively detecting 
cardiac phenomena of claim 38 wherein said means for 
generating bispectral values comprises a means for 
computing the Fourier transform of the third order 
crosscorrelation function of said filtered signals. 
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STATEMENT UNDER ARTICLE 19 



Claims 1, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 
38, 40, 45, 48, 49, 50, 51, and 52 of the above-referenced 
application have been amended to more clearly describe the 
subject of the invention and to make them consistent with 
the amended claims in the priority bnited States 
application (which claims have been amended to cure 
indef initeness problems), and to correct obvious errors. 

Claims 53, 54, 55, and 56 have been added to clearly 
protect various means for generating bispectral values in 
applicant's system. 

Replacement sheets (pages 39-42 and 46-49) are 
enclosed for pages containing amended, or, new claims. 
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